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Abstract 

We study the frictional drag due to Coulomb and phonon mediated electron-electron interaction 
in a double layer electron system exposed to a perpendicular magnetic field. Within the random 
phase approximation we calculate the dispersion relation of the intra Landau level magnetoplas- 
mons at finite temperatures and distinguish their contribution to the magnetodrag. We calculate 
the transresistivity pDrag as a function of magnetic field B, temperature T, and interlayer spacing 
A for a matched electron density. For A = 200 nm we find that PDrag is solely due to phonon 
exchange and shows no double-peak structure as a function of B. For A = 30 nm, porag shows the 
double-peak structure and is mainly due to Coulomb interaction. The value of PDrag is about 0.3 
f2 at T = 2 K and for the half-filled second Lanadau level, which is about 13 times larger than the 
value for A = 200 nm. At lower edge of the temperature interval from 0.1 to 8 K, pDrag/T 2 remains 
finite for A = 30 nm while it tends to zero for A = 200 nm. Near the upper edge of this interval, 
PDrag for A = 30 nm is approximately linear in T while for A = 200 nm it decreases slowly in T. 
Therefore, the peak of PDrag/T 2 is very sharp for A = 200 nm. This strikingly different magnetic 
field and temperature dependence of PDrag ascribe we mainly to the weak screening effect at large 
interlayer separations. 
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I. INTRODUCTION 



Frictional drag between two spatially separated two-dimensional (2D) electron systems 
is a powerful tool to probe electron-electron (e — e) interaction and is of current intensive 

nQ 

research interest both experimentally and theoretically |H,|2j. The drag effect manifests itself 
when a current with a density J\ driven along the layer 1 induces, via momentum transfer, 
an electric field E 2 in the layer 2 under the condition that the layer 2 is an open circuit 
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Much experimental work was done in the magnetic field-free case to measure the drag 
effect in a coupled 2D and 3D electron system j^, between two 2D electron systems (2DES) 



two 2D hole systems (2DHS) 13J, Il4j . and also in a coupled 2DES 



2DHS 13j, [15|. Theoretical work on the frictional drag has 



3een devoted to explain different 



momentum transfer mechanisms in a coup 



2DES-2DES 
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7J, 2DES-2DHS 
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33J, and recently in a 2DHS-2DHS 



Also, plasmon-enhancement of the dra g du e to many-body correlations at high tem- 



peratures has been studied theoretically [28, 30] and experimentally 



Two possible 



interlayer e — e interaction mechanisms have been considered: One is the direct electro- 
static Coulomb scattering and the other is the phonon-mediated effective interaction. The 
Coulomb drag calculations predict a T 2 temperature dependence for the transresistivity 



p Dra and a strong inter- layer spacing (A) dependence as A -4 [23L |25|. However, the ex- 
perimental results do not support the pure Coulomb mechanism [jl lsl \v\ . The observed 
drag has a peak when the electron densities match, n\ = n 2 , and the scaled transresistivity, 
p Drag /T 2 , deviates from a constant value demonstrating a peak as a function of T well below 
Bloch-Griineisen temperature. Also, the measured p Drag for separations at least up to 50 
nm exhibits almost no dependence on A and its value beyond A = 50 nm is much bigger to 
be accounted for by Coulomb interaction alone. This led Gramila et al. |7[ to propose an 
acoustical phonon-mediated scattering mechanism for the drag. This mechanism was studied 
extensively in Refs. The theoretical results are satisfactory to explain the main 

experimental findings. The phonon-mediated drag due to the bare electron-phonon (e — p) 
coupling diverges when the transferred energy, hu, from the layer 1 to the layer 2 is close to 
hsq (if is the transferred in-plane momentum). This divergence is eliminated by taking into 
account other scattering mechanisms such as screening and/or the finite phonon mean free 
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path (32I. It is found that, despite the weak e — p coupling, the effective phonon-mediated 
e — e interaction remains strong and competes with the Coulomb interaction. In general, 
e — p scattering in 2DESs is qualitatively different between in the small angle scattering 
(SAS) temperature region, T <C tiskp, and in the large angle scattering (LAS) temperature 
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region, T ~ hskp, where s is the speed of sound and kp the Fermi wave vector 
the SAS region, %uj < T and this results in the power law dependence of porag oc T 6 . On the 
other hand, in the LAS region, tkv ~ 2hskp and pDrag varies linearly in T. As a consequence, 
PDrag/T 2 has a peak at the crossover from the SAS to LAS temperature region at the peak 
temperature T pea k < hskp- In the LAS temperature region the momentum transfer is most 
efficient when momenta of the emitted and absorbed phonons are determined by the same 
Fermi wave vector k\p = k 2 F, this results in the peaked transresistivity when n\ = n 2 . In 
the SAS region there is no cutoff related to kp that pnrag shows a monotonous behavior as 
a function of n 2 . 

The frictional drag has been also investigated in the presence of a perpendicular 
magnetic field. The drag rate has been measured experimentally in a coupled 2DES 



j. Also, the response of 



2DHS 13] as well in a 2DES-2DHS 
composite fermions to large wave vector scattering has been studied through phonon mag- 
netodrag measurements |44|. In addition, recent measurements reveal a new regime of the 



magnetodrag in a coupled 2DES at mismatched densities n\ 7^ n 2 in which the polarity of 
the drag voltage is opposite that normally found for electron systems 3- Theoretical 
inve stig ations were mostly carried out with incorporating the direct Coulomb mechanism 
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471 . 148L |4j| . These calculations show reasonable agreement with the experiments 27J, |3j| 
that at low temperatures porag demonstrates an analog of the Shubnikov-De Haas oscilla- 
tions with amplitudes by two order of magnitude enhanced compared to the zero field drag 
signal. At the matched densities the measured magnetodrag has a double-peak structure 
in magnetic field in the inter-quantum Hall plateau regions [3, Q|- This observation is 
in agreement with the earlier calculations of the Coulomb magnetodrag by B0nsager et al. 
who predicted this behavior caused by the interplay of screening and Landau quan- 



;ization. However, the calculations by Wu et al. |48[ and the experiments by Hill et al. 
371 1 show oscillations without the double peak structure. The critical test of the theory was 

n 

performed in a recent experiment by Lok et al. |45|], which confirms that the transresistiv- 
ity does not show the predicted double peak structure for the spin split Landau levels and 
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the double-peak structure at higher filling factors is not caused by the screening effect |43j. 
Notice also, as against to the zero magnetic field case, the scaled transresistivity pDrag/T 2 
due to the pure Coulomb mechanism shows a peaked temperature dependence in the finite 
magnetic fields 47 1. In addition, the Coulomb magnetodrag calculations j3, Isoj l 
reveal new features associated with interlayer magnetoplasmons, that was observed also 
experimentally P|. 

As we described in the previous paragraph, the phonon-mediated drag is the dominant 
mechanism in the zero magnetic field case in a coupled 2DES, except for very closely spaced 
layers [3,0, [ill Further, it was reported recently that the effective phonon and 

the Coulomb contributions to drag were comparable in a coupled 2DHS even with closely 
spaced layers. The phonon-mediated drag effect has been also investigated experimentally 
in the finite magnetic fields 13, [3] where it was found also that the phonon contribution 
was prevalent in the coupled 2DES, 2DHS, and 2DES-2DHS with large barriers. Up to 
now, however, it is rare to find theoretical works on the frictional drag mediated by phonon 
exchange in a finite magnetic field. The exceptions are recent treatments of the phonon 
dra g at the Landau level filling factor v = 1/2 by Chern-Simons composite fermion theory 

The aim of the present paper is to investigate the direct Coulomb and the effective 
phonon-mediated frictional drag in a coupled 2DES exposed to a perpendicular magnetic 
field theoretically. For this system we calculate the dispersion relation of the intra-Landau 
level magnetoplasmons by taking into account finite temperature and distinguish the mag- 
netoplasmon contribution to the magnetodrag. We calculate porag as a function of magnetic 
field B, temperature T, and interlayer spacing A, and determine the relative contributions 
of Coulomb versus phonon-mediated e — e interaction to magnetodrag 

This paper is organized as follows. In the next Sec. II we describe briefly the theoretical 



32, 



47. In 



method we use to calculate the transresistivity, which follows mainly to Refs. 
Sec. Ill we present and discuss our calculations. Finally, the results are summarized and 
the conclusions are given in Sec. IV. 
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II. THEORY 



Our theoretical model consists of two parallel electron layers exposed to a perpendicular 
magnetic field in the z-direction, separated by a distance A from center to center, each having 
the electron densities n\ and n 2 and the layer extensions d\ and The experimentally 
measured quantity is the transresistivity defined as porag = E 2 / J\. The main theoretical 
framework to calculate porag is the Fermi liquid theory where the interlayer e — e interaction 
is treated perturbatiyely. The expression for porag can be derived by using the linearized 
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tzmann equation |25( or the memory function formalisms starting from the Kubo formula 
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32j. Direct calculations show that in our model the transresistivity is represented 
in the standard form 

h 2 1 ^ 2 r°° duo \W 12 (q,uj)\ 2 
pDrag ~ 2eV 2 ri| 9 Jo 2vr \e(q,u;)\ 2 

smh. 2 (hw/2T) 1 ' 

where A is the normalization area, Wij(q, u) = J2r ^ / iT(9) u ) ^ ne total unscreened interlayer 
interaction matrix element = 1,2 are the layer indices), e(q,u) the screening function, 
and Xi{Qi u ) the irreducible intralayer electron polarization functions (the interlayer electron 
polarization is neglected). The index T refers to the type of e — e interaction. We consider 
direct Coulomb (T = C) e — e interaction and effective e — e interaction, mediated by the 
exchange of piezoelectric (T = PA) and deformation (T = DA) acoustical phonons as well 
as the exchange of polar optical phonons (T = PO). Phonon mediated e — e interaction 
appears in second order perturbation theory with respect to the bare e — p coupling. The 
unscreened interaction matrix element of this process is visualized by the diagram in Fig. 
The free phonon propagator 



D r (Q,(j) = 2^/r 1 



r\ 2 



(2) 



with Q = (~c?,q z ), Uq ,da = sQ and Uq° = <jJ PCn and a summation over phonon momenta 
corresponds to the internal phonon line (the dashed line in the diagram). Tq is the lifetime 
of the phonons. The vertex part of this diagram (the solid dots in Fig. corresponds to 



r/ = Jb-twm (3) 
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FIG. 1: Diagram for the effective phonon mediated e — e interaction. Here Ei, k x , and a are the 
electron energy, the momentum in the Landau gauge, and the spin, respectively. 

where B T (Q) are the bare e — p coupling functions j3] and = J dzpi(z)e^ z ^ d are 

the form factors in z-direction (the remaining in-plane part of the form factor enters the 
definition of u;)). We assume that electron scattering takes place in the lowest electron 
subbands of infinitely high quantum wells with widths d\ = d,2 = d and the subband density 
functions P2{z) = Pi(z + A) in both layers do not depend on the subband index. In GaAs 
the PA interaction matrix element is given by the effective piezoelectric modulus which 
depends only on the polarization of the phonon and its direction of propagation. Since in 
this description of frictional drag we are interested in the average relaxation characteristics 
of the system, the anisotropy of PA interaction can be ignored. We describe PA interaction 
within the scope of the isotropic model [3] that leads to a scalar PA interaction constant. 
We assume also that all elastic parameters of the sample are the same and phonons are not 
reflected from the interfaces separating different materials (for instance between GaAs and 
AlGaAs). 

Taking the summation over q z and including also the Coulomb contribution, we obtain 
for the total unscreened interlayer e — e interaction matrix elements 

K q { Tpo uj po — LO 

+— tiAad) 

t pa a 

+^ Us* + TT-^l F ^ ad) ) ) (4) 

tda V 2 (sp PO ) a J J 



U Q i 



where n = 13.1 |54( is the GaAs static dielectric constant, hp PO = ^2mhtu PO , u PO the 
polar optical phonon frequency, and a = Jq 2 — (7 + 2s7 r ) 2 - m ^ ne li m it of infinite phonon 



lifetime a = Jq\ — uj 2 / s 2 if q\ > uj 2 / s 2 and a = —iyQ± — uj 2 / s 2 if q\ < uj 2 / s 2 with the 
positive square root branch. The form factors F are given by 

[ F 2 (0,ifi = j, 

F 2 (£) = J dz 1 dz 2 p(z 1 )p(z 2 )e^-^ 2l ,F 3 = dj dzp\z). (5) 

The following nominal scattering times are introduced in Eq. (0J) 

l_ 2upo J_ _ B? A p PO PA _ h{e(5f 

r c p PO a t pa 27m s 2gs 

1 1 B DA r? h" 2 

^- = 2a F u PO , J- = 2L*«L, = ( 6 ) 
r PO r DA vrft s 2gs 

where olf is the Frolich interaction constant, 5 and e/3 the deformation and piezoelectric 

phonon potential constants, g the crystal mass density, and a* the effective Bohr radius. For 

GaAs we have T q ~ 0.024 ps and the numerical values Tpo ~ 0.14 ps, T PA ~ 8 ps, T/x4 ~ 4 



ps take we from 
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We are interested in the experimental situations when the Landau levels are fully resolved, 
accordingly the scale of u over which a 2DES responds to an external perturbation is given 
by the width of the Landau band. Therefore, when u <C uipo, the optical phonon propagator 
D PO (q , uj) is never on the mass surface, and only virtual optical phonons can contribute to 
the magnetodrag in this regime. The PO phonon mediated contribution to total interaction, 
(q, uj) (the second term in Eq. (Jlj)). at uj = 0, behaves like the Coulomb contribution, 
W£j(q,uj) (the first term in Eq. And, because of the small t c /t po coupling, simply 

results in about 17% renormalization of W^(q,u). From Eqs. and (JHJ) one can see that 
W^(q,uj) is large when qA < 1. In drag experiments A is larger than d so we have usually 
qd <C 1. The PA and DA phonon contributions, W^ A ' DA (q,uj) (the third and forth terms in 
Eq- ©), behave like a delta function because of the small phonon couplings t c /t pa ps 3- 10~ 3 
and t c /t pa 1.5 • 10 -3 : W? A ' DA are negligible with respect to the Coulomb term W?(q, uj) 
for all over uj and q except the close neighborhood of uj = sq where a is very small and we 
have again ad <C 1. Thus, the zero-limit approximation of the form factors F is well justified 
for both Coulomb and phonon mediated e — e interaction and we take Fi t2 (^) ~ -^1,2(0) in 
our numerical calculations. 

In the random phase approximation the dielectric tensor is obtained from the solution of 
a matrix Dyson equation for the dynamically screened interlayer interaction [55] and has the 



form Eij(q,u) = 8ij — Wij(q,uj)xj(q,u) (assuming no summation by the repeating indices). 
The screening function e(q,u) is the determinant of the dielectric tensor and is given by 

e(q,u) = (1-W n (q,u)xi(q,u)))(l-W 22 (q,uj)xi(q,u)) 

-W? 2 (q,uj)xi(q,uj)x2(q^) (7) 

In this approximation yfq. uo) is determined by the bubble diagrams with the exact vertex 
part. B0nsager et al. [47] have shown that when the Landau levels are clearly resolved, 
corrections to the vertex part are small. In the first approximation they can be neglected 
and xio.) u ) is given by 



x 



JmG?{E) [pf (E + hu) + G${E - hu)) , (8) 



Ib is the magnetic length and fr the Fermi distribution function determined by the chemical 
potential Ep. The bare vertex functions Qu> are given by the gauge invariant part of the 
in-plane form factor 

Q w it) = (/!//'!) 1/2 e-W-'Wf 1 ' (t) , for I < I' (9) 

where L\t 1 (t) is the associated Laguerre polynomial. For V > I one can exploit the symmetry 
relation Q w (t) = {-l) 1 - 1 ' Qw (t) • In Eq. © G?' A (E) are the electron retarded and advanced 
Green functions dressed by electron-impurity (e — i) interaction. We treat e — i scattering 



in the short-range impurity model 
Uemura 



and use the Green functions obtained by Ando and 



571 ] in the self-consistent Born approximation 
Gf A {E) = 



E — Ei + J (E — Ei) 2 — Tq 



(10) 

In this approximation the half-width To = \J^Tiujb^: of the Landau level E\ — (I + l/2)huB 
is independent of the Landau index I {ub is the cyclotron frequency and r is the transport 
relaxation time that determines the current in the drive layer via the mobility \i = er/m*, 
m* is the electron effective mass). The chemical potential Ep(n, B, T, /x) in distribution 
functions is determined implicitly from the electron density 



n 



Here we use the Green functions obtained by Gerhardts 5^| in the improved self-consistent 
Born approximation 

^ _ / (E-E t ) 2 \ 



G t (E) = -i^— exp -2 

in V 1 



X 



o \ 1 o 

E-Ei 



1 + Er, 



n^^jj (i2 » 



which correspond to the Gaussian density of states without the unphysical edges of Landau 
bands. Erf(z) is the error function. 

Thus, Eqs. (JTJ) and (P IT^) determine the transresistivity of the double layer 2DES exposed 
to a perpendicular magnetic field and coupled via direct Coulomb and effective phonon 
mediated e — e interaction. In this treatment the Landau levels acquire broadening due 
to e — i scattering and each broadened state has its weighted contribution to drag. The 
electron energy remains dispersionless and the distribution functions do not depend on the 
position of the Landau oscillator center. For this reason porag does not depend on the gauge 
non-invariant electron momenta k, as it should be, but also after performing the summation 
over all k, no interference occurs in Eq. between different Fourier components "Tf both 
of the Coulomb 59] and the e — p interaction potentials. It should be noticed also that only 
such non-trivial treatment of e — e, e—p, and e — i scattering allows the existence of the intra- 
Landau band magnetoplasmons and e — p scattering processes with a finite energy transfer 
hu 0. This is especially important in the regime of T, To -C %ujb under the consideration 
here. Without impurity broadening of the Landau levels, the phonon mediated magnetodrag 
in this regime is only due to exchange of pure virtual phonons with Tiuo = and q ^ 0. But 
this contribution is small with respect to the Coulomb magnetodrag for all phonon modes. 
In contrast to this, for finite hu, acoustical phonon mediated e — e interaction diverges in 
uj = sq. This divergence enhances strongly the phonon contribution to magnetodrag and 
makes it dominant at large interlayer separations. 



III. RESULT AND DISCUSSION 



Below we present our calculations of pDrag as a function of B, T, and A in a symmetric 
double layer quantum well system with a matched electron density n = 2.5 • 10 15 m -2 . We 
neglect the spin splitting and assume that the Landau levels are fully resolved and the odd 
filling factors correspond to a half-filled Landau levels. The spin degeneracy results in an 
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0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 2: Intra Landau level magnetoplasmon dispersion relation for the adopted parameters ob- 
tained from the zeros of the real part of the screening function. Dotted line shows zeros of the real 
part of the screening function for a single layer. 

additional factor of 4 in Eq. {I}. It is known that at temperatures below Bloch-Griineisen 
temperature the frictional drag is dominated by the PA phonons in the case of B = 
j&J, and the PO contribution is weak as we describe in the previous Section. Here we 
consider mainly PA phonon mediated e — e interaction. The inclusion of calculations for 
DA phonon mediated e — e interaction is straightforward. We take into account only the 
intra-Landau level magnetoplasmon modes and, in the low temperature regime (T <C ujb) 
under consideration here, we neglect the interlevel magnetoplasmon contributions to the 
magnetodrag. The intra-Landau level collective excitations £0j exist when a Landau level is 
partially filled. In Fig. Elwe plot the intra-Landau level magnetoplasmon dispersion relation 
obtained from the zeros of the real part of the screening function, Re e(q,u(q)) = 0, 
with the imaginary part describing the damping. The spectra have two collective modes 
with the finite values (0) and uj + (0) close to r /2 and 3r /2 in the limit of q — > 0. 
When B varies in the second inter-plateau region from 3 to 4.5 T for T = 2 K or T 
varies from 1 to 8 K for B = 4 T, oj~ (0) (u + (0)) is changed slightly, within 1% (5%). 
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The values of uj~ (0) and u + (0) increase monotonically in T in this temperature interval. 
In this magnetic field interval, they approach closely each other at B corresponding to 
the half-filled Landau level. The whole spectra are located in a finite range of q close to 
zero so that the magnetoplasmons have momenta q£s < 1 for the choice of parameters in 
Fig. |21 For such momenta the phonon energy hsq (the dashed-dotted line in Fig. |2J) lies 
below the magnetoplasmon energy near the lower edge of Landau band. Therefore, the 
phonon mediated e — e interaction effect on the magnetoplasmon spectrum is negligible. 
For energies Too well below 2r (r ~ 6.5 K, Tq/Tiujb ~ 0.08 for B = 4 T), both the 
splitting and broadening of each magnetoplasmon mode of the spectra are determined by 
the strong damping Im x{Qi ^(q))- The poles corresponding to these collective modes are 
moved into the complex tu-plane and are located far from the real axis. Therefore, the two 
modes of the spectra are strongly mixed and damped and cannot be interpreted as the 
normal modes of the system. We note that in this strongly dispersive system, the group 
velocity v g = du(q)/dq shows anomalous behavior: It becomes infinite and changes its sign 
at certain points. This should not cause an alarm because in the regions of anomalous 
dispersion with Re s(q,u(q)) = 0, the approximation Au(q) « v g q is not valid and the 
group velocity is not a useful concept [6lJ . Noticed also that the situation when the group 
velocity changes its sign in a finite wavenumber arises also, for instance, in the single particle 

and 2D 63, Q] electron systems exposed to a 



spectra of magnetic edge states in ID 
nonhomogeneous magnetic field. Towards the upper edge of the Landau band the damping 
decreases and becomes zero for Hu > 2Tq. Thus, in this model of the Landau band with 
sharp unphysical edges, the magnetoplasmons with energies Tiuj > 2To have no contribution 
to the magnetodrag. Near the upper edge of Landau band, the spectra are splitted into 
the symmetric and antisymmetric magnetoplasmon modes and this is due to the interlayer 
coupling, determined by the spacing A. One can see that for A = 30 and 200 nm the 
magnetoplasmon dispersion shows different behavior in the close vicinity of hu = 2Tq. In 
this region damping is vanishingly small, and the spectra describes long lifetime bound states 
of magnetoplasmons in a double layer system with a binding energy given by the splitting 
of the magnetoplasmon modes. For A = 30 nm the binding energy is about 1 K near the 
upper edge of the Landau band, and it decreases strongly with A. For A = 200 nm Coulomb 
interlayer e — e interaction is obviously very weak, the binding energy is negligibly small, and 
the magnetoplasmon dispersion becomes degenerated. We note that the spectrum of inter 
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FIG. 3: Transresistivity as a function of the interlayer separation, normalized to the value of total 
transresistivity at A = 20 nm. The contributions of the pure Coulomb and the pure phonon 
mediated couplings are shown together with their combined total contribution. 

Landau level magnetoplasmons near the upper edge of the Landau band has been calculated 



50, where 



analytically for pure Coulomb e — e interaction by Khaetskii and Nazarov in Ref. 
the analogous dispersion curves for higher cyclotron harmonics was obtained. 

In Fig. El we plot the transresistivity as a function of the interlayer spacing A for pure 
Coulomb e — e interaction, W^(q,u), and for pure PA phonon mediated e — e interaction, 
WE (q, uj). In addition, we plot the total porag due to their combined effect, Wij(q,u) = 
W£j(q,ui) + W? A (q,ui). In all three cases we calculate the screening function by taking 
Wij(q,cj) = W^(q,uj) + W[j (q,ui). It is seen that the Coulomb mechanism dominates for 
A < 30 nm with an approximately exponential decrease of porag m A while the phonon 
mechanism dominates for A > 80 nm with an approximately logarithmic decrease. They 
have comparable contributions to magnetodrag for the intermediate A, showing the equal 
contributions at about A = 50 nm. 

In Fig. 0] we plot the magnetic field dependence of total transresistivity due to Coulomb 
and PA phonon mediated e — e interaction. The transresistivity shows Shubnikov-De Haas 
oscillations as a function of B both for the A = 30 and 200 nm spacings. For A = 200 nm we 
find that PDrag{B) has perfectly symmetric oscillations with strongly increasing amplitudes 
in B. We ascribe this to weak screening at large interlayer separations. For A = 30 nm, 
in agreement with the calculations by B0nsager et al. |4f| for pure Coulomb interaction, 
PDrag(B) demonstrates a double peak structure around the filling factor v = 3, caused by 
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FIG. 4: Total transresistivity as a function of magnetic field plotted for the A = 30 and 200 nm 
interlayer separations. The data for A = 200 nm has been multiplied by a factor of 26. 



strong screening. At low fields screening is relatively weak and, except for some asymmetry, 
no dip is visible for v > 4. However, screening has a flattening effect on the oscillation 
amplitudes in different inter-plateau regions for A = 30 nm. The phonon contribution to 
the total magnetodrag for A = 30 nm is about 5 to 10% in the second inter-plateau region, 
depending on B. For A = 200 nm the Coulomb contribution to PDrag{B) is negligible. It is 
seen from Fig. HJthat near v = 3 the total porag is about 0.3 Q for T = 2 K and A = 30 nm, 
and this is approximately 13 times larger than pm-ag for A = 200 nm. The magnitude of the 
calculated pDrag and its strong decrease with increasing A are in good quantitative agreement 
with the experimental findings by Rubel et al. In this experiment, however, dips are 
also observed in the middle of Landau bands for samples with large interlayer separations. 
The lack of such dips in our calculations for A = 200 nm might be an additional indication 
that these dips in the experiments should be attributed rather to the spin than the screening 
effect BQ, 3- 

In Fig. Elwe present the temperature dependence of porag and pDrag/T 2 for two different 
values of B and for A = 30 and 200 nm in the temperature range from 0.1 up to 8 K. 
Although porag/T 2 shows a peaked temperature dependence for both A = 30 and 200 nm, 
PDrag demonstrates completely different behavior for A = 30 and 200 nm as function of 
T. For A = 30 nm the magnetodrag is mainly due to the Coulomb mechanism and is 
in agreement with the calculations by B0nsager et al. j3| for pure Coulomb interaction. 
At lower temperatures pDrag has a subquadratic temperature dependence and porag/T 2 
remains finite at low T. The transferred energy fiu is restricted by temperature, Tiuo < T, 
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FIG. 5: Total transresistivity (the left axis) and scaled transresistivity (the right axis) versus 
temperature for B = 3 and 4 T and for A = 30 nm (thin lines) and 200 nm (thick lines). The data 
for A = 200 nm has been multiplied by a factor of 26. The electron density is n = 2.5 • 10 15 m~ 2 , 
the phonon mean free path L/j = 100 //m, and the mobility fi = 25 m~ 2 V -1 s _1 . 



and the transferred momentum is, independently of hu, restricted by the interlayer spacing, 
q < A -1 (in quantizing magnetic fields A -1 is always smaller than i]} 1 ). This is clearly 
seen in Fig. El where we plot the integrand over hu in Eq. (JTJ for A = 30 nm after taking 
the sum over if. Even for relatively large T = 1 K the integrand is strongly decreasing 
function of to. Therefore at low temperatures one can use the small uj and the small q 
limit of the polarization function ^(?,w). In this limit both Im xiQ^) an d Im e(q,u) are 
linearly vanishing functions with uj while Re e(q,u) has a large maximum, determined by 
the electronic compressibility. By replacing the upper limit of the integration over Two in 
Eq. (JTJ) with T and taking the static limit of the screening function out of the integration, 
one can obtain (cf. Fig. |HJ) 



p Drag OC T 2 / 
J 



T dhuj (Imx (M IT) 2 

OC / 

T sinh 2 (^/2T) 
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FIG. 6: Integrand over uj in Eq. ((J) after taking the sum over ~q for A = 30 nm and for different 
temperatures. The data for T = 1 K has been multiplied by a factor of 40. 

At high temperatures all %uo < 2r contribute to the integral in Eq. (JIJ). One can see from 
Fig- El that for T = 3, 6, and 8 K the integrand decreases very slowly in the region of small 
uj. However, at energies close to Tiuj~ (0), when the first magnetoplasmon mode is developed 
(cf. Fig. |2J), the integrand starts to increase slowly. Further, at energies close to huo + (0), 
when the second magnetoplasmon mode is developed (cf . Fig. EJ , the integrand starts to 
increase strongly. The integrand reaches its well pronounced maximum, with an increasing 
height in T, just before it drops to zero at the upper edge of the Landau band Tiu = 2T . 
This behavior of the integrand results in the peaked temperature dependence of porag/T 2 
shown in Fig. El with an approximately linear temperature dependence of porag oc T at high 
temperatures, observed experimentally by Hill et al. jjj. 

The situation is completely changed for A = 200 nm where the magnetodrag is mainly 
due to phonon exchange. At low temperatures porag has stronger dependence on T than 
it has at the small interlayer separation. At high temperatures porag decreases slowly with 
increasing T. This results in the very sharp peak of PDrag/T 2 as a function of T for A = 200 
nm (Fig. EJ). This strikingly different behavior of pDrag as a function of T for different A is 
a consequence of the fact that the main contribution to the phonon mediated drag comes 
from q uj / s for which effective unscreened e — e interaction diverges in the limit of infinite 
phonon lifetime. Now the integrand over u in Eq. (0) is proportional to gi(t) = t (Li(t)) 2 e~ l 
with t = (oj/s) 2 £ 2 B /2 and in the whole temperature range the main contribution to phonon 
magnetodrag makes t = t\ ~ 1 (t\ are the zeros of the derivative of gi{t)), i.e. the finite 
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hu ~ hsq ~ hs/is ~ To- The small u and the small q approximation is not valid for e(q, u) 
and it is rapidly decreasing function with u and the screening effect is weak in the phonon 
magnetodrag (cf. the curves corresponding to A = 30 and 200 nm in Fig. HJ). We obtain 
the following formula as an approximation to describe the temperature dependence of the 
transresistivity in this regime 

\- 1 U lWujt^y^ 2 {lm X (ti,yi)) 2 
PBrt-lfM \ £ {t hVl )\* sinh 2 (^/2^) U4j 

where d = T/Tq and y = Tioj/T$. For B = 4 T we have v = 2 and the main contribution 
to the sum makes the outermost filled Landau level with I = 1 and ti = 2 + y/3. In GaAs 
2ms 2 ~ 0.23 K and r ~ 6.5 K for p = 25 m 2 V _1 s _1 , so we have t ~ (1.5y) 2 . This gives 
Tuo\/2 « 4.2 K as a crossover temperature of the temperature dependence of pDragiT). At 
temperatures below Uui/2, sinh 2 (frci;i/2T) increases strongly with T and this determines 
mainly the behavior of PDragiT). It is seen from Fig. Elthat for A = 200 nm porag/T 2 tends 
to zero at low T in contrast to its finite value for A = 30 nm. At temperatures above huui/2, 
sinh 2 (/iu;i/2T) oc T 2 , and the screening function e(qi,coi) for such large arguments is small 
and increasing in T. This behavior together with the T _1 prefactor in Eq. (0) determines the 
temperature dependence of purag at high T. One can see, however, from Fig.Elthat the actual 
peak position of pDrag{T) is shifted to the left from huj\/2. For A = 200 nm the integrand 
over io in Eq. Q after taking the sum over if manifests itself to three strongly pronounced 
peaks near y = 0.45, 1.25, and 2, and their positions are relatively stable with respect to 
temperature (see Fig. |7|). Notice, as against the situation of the Coulomb magnetodrag, 
even for small T = 1 K the integrand is approximately zero in the region of small uj and has 
strongly pronounced peak at large u. The main highest peak of the integrand near y = 1.25, 
as we already discussed, is related to the maximum of gi(t) at tf = 2 + v^3 corresponding 
to y = yi = 1.29. The left small peak of the integrand near y = 0.45 is related to another 
maximum of the function gi(t) at t\ = 2 — \/3 (y{ = 0.35) and also to the maximum of 
the function ^o(^) at to = 1 (Vo — 0.67). This second peak contribution to PDragiT) shifts 
the actual crossover temperature of pDragiT) to the left so that T pea k ~ 3.6 K (Fig. EJ). 
One can see, however, that the positions of the first two peaks of the integrand is not 
far from the magnetoplasmon energies hu^ (0) . Therefore the magnetoplasmons contribute 
also to the formation of these peaks to some extend. Nevertheless, it is clear from the very 
weak dependence of the total porag on A (Fig. |3J) that the magnetoplasmon contribution 
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FIG. 7: Integrand over uj in Eq. (^Q) after taking the sum over ~q for A = 200 nm and for different 
temperatures. 

to porag for A = 200 nm is a small correction to the phonon contribution. The third 
peak of the integrand in the very close vicinity of the upper edge of the Landau band is 
a pure magnetoplasmon effect (cf. Fig. Eland Fig. [7J). Immediately below the upper edge 
of Landau band the magnetoplasmon damping is very small and this results in the very 
sharp peak of the integrand with a relatively small contribution to PDrag{T). Thus, at 
high temperatures we obtain that for A = 200 nm porag decreases slowly with T (Fig. [5J). 
This extraordinary temperature dependence of porag is mainly due to the weak screening 
effect at large interlayer separations and the screening is weak because the phonons with 
large energies and momenta realize e — e interaction between remote layers. Notice that 
the T~ l dependence of PDrag{T) for the pure Coulomb interlayer interaction in one of the 
typical scattering regions was also reported by Khaetskii and Nazarov |50j]. As far as we 
know at present no magnetodrag measurements are available in the above regime and the 
experimental test of the temperature dependence of phonon magnetodrag is critical. 



IV. SUMMARY AND CONCLUSIONS 



We have calculated the transresistivity between spatially separated electron layers in a 
perpendicular magnetic fields. We take into account both direct electrostatic Coulomb and 
effective phonon mediated e — e interaction. For this system we calculate the dispersion re- 
lation of the intra-Landau level magnetoplasmons within the random phase approximation 
at the finite temperatures and distinguish the magnetoplasmon contribution to the magne- 



17 



todrag. We find the strikingly different magnetic field and temperature dependence of the 
transresistivity for the small and large interlayer separations. When A = 30 nm, the transre- 
sistivity shows a slight dip as a function of B in the middle of the second spin-degenerated 
Landau band with the total pDrag about 0.3 Q (the phonon mean free path is 100 /im, the 
mobility p. = 25 m 2 V -1 s _1 , T = 2 K, and n = 2.5 -10 15 m -2 ), which, depending on B, 
includes about 5 to 10% of the phonon contribution due to piezoelectric interaction. When 
A = 200 nm, the Coulomb contribution to magnetodrag is negligible. And, porag shows 
no dip and the total drag around filling factor v = 3 is about 13 times less than the value 

or A = 30 nm. This is in good agreement with the experimental findings by Rubel et al. 

3^| . The scaled transresistivity porag/T 2 has a peaked temperature dependence for both 
A = 30 and 200 nm. At low T, porag/T 2 remains finite for A = 30 nm, while for A = 200 
nm it tends to zero. At high T, porag is approximately linear in T for A = 30 nm while for 
A = 200 nm it decreases slowly with T. Therefore, the peak of PDrag/T 2 as function of T 
is very sharp for the large separation A = 200 nm. We ascribe this to the weak screening 
effect at large interlayer separations where the phonon mechanism dominates with the large 
transferred energy and momentum. 
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